Collective dynamics of liquid aluminum probed by 
Inelastic X-ray Scattering. 
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An inelastic X-ray scattering experiment has been performed in liquid aluminum with the purpose 
of studying the collective excitations at wavevectors below the first sharp diffraction peak. The high 
instrumental resolution (up to f .5 meV) allows an accurate investigation of the dynamical processes 
in this liquid metal on the basis of a generalized hydrodynamics framework. The outcoming results 
confirm the presence of a viscosity relaxation scenario ruled by a two timescale mechanism, as 
recently found in liquid lithium. 

PACS numbers; 6f.lO.Eq, 67.40.Fd, 67.55. Jd, 63.50.+X 



I. INTRODUCTION 

In the last two decades, the dynamical properties of 
liquid metals have been widely investigated aiming at 
the comprehension of the role of the mechanisms under- 
lying the atomic motions at the microscopic level. In par- 
ticular, in the special case of alkali metals, it is known 
that well-defined oscillatory modes persist well outside 
the strict hydrodynamic region, down to wavelengths 
of a few inter-particle distances, making these systems 
an ideal workbench to test the different theoretical ap- 
proaches developed so far for the microdynamics of the 
liquid state. By the experimental point of view, it is 
worth mentioning the pioneering inelastic neutron scat- 
tering (INS) experiment by Copley and Rowe ^ in liquid 
rubidium, while more recently a lot of experimental ef- 
forts in performing more and more accurate experiments 
have been done: INS investigations have been devoted 
to liquid cesium [Q, sodium hthium Q, potassium 
lH and again rubidium [|j . Many numerical studies have 
been reported on the same systems, allowing to overcome 
experimental restriction such as the {Q — E) accessible 
region when dealing with the collective properties (with 
the density fluctuation spectra), and the simultaneous 
presence of coherent and incoherent contribution in the 
INS signal, the paramount experimental technique in this 
field up to a few years ago. Moreover, numerical tech- 
niques allow to access a wider set of correlation function 
while the inelastic scattering experiment basically probe 
the density autocorrelation function. 

On the theoretical side, tanks to the development of 
tools such as the memory function formalism, the relax- 
ation concept, the kinetic theory, a framework has 
been established in order to describe the behavior of the 
afore mentioned correlation functions. In this respect, 
among the most important task, there is the idea that 
the decay of the density autocorrelation function occurs 
over different time scales. In particular a first mecha- 



nism is related to the coupling of the density and the 
temperature modes (thermal relaxation), while a second 
process involves the stress correlation function (viscosity 
relaxation) . The latter mechanism proceeds through two 
different relaxation channels, active over different time 
scales: a first rapid decay, that has been ascribed to 
cage effects due to microscopic interaction of each atom 
and the cage of its neighboring atoms, followed by a 
slower process related to the long time rearrangements 
that drives the glass transition in those systems capable 
of supercooling. This structural relaxation process has 
been widely studied in the mode coupling formalism and 
recently remarkable efforts have been developed to set 
up self consistent approaches |^ . This theoretical frame- 
work has been tested on a number of numerical studies 
(see for example [0 and references therein). Indeed, 
the lineshape of the density correlation spectra extracted 
by INS (in those systems where it is allowed in the sig- 
nificant (Q — E) region) did not allow to discriminate 
between different models including all the relevant relax- 
ation processes at the microscopic level |^ . 

In the recent past, the development of new synchrotron 
radiation facilities opened the possibility of using X-rays 
to measure the S(Q, u) (which is proportional to the scat- 
tered intensity) in the non-hydrodynamic region; in this 
case the photon speed is obviously much larger than the 
excitations velocity, and no kinematic restriction occurs. 
Moreover, in a monatomic system, the Inelastic X-ray 
Scattering (IXS) cross section is purely coherent and so 
it is directly associated with the dynamic structure fac- 
tor. 

Recently, an accurate IXS study on liquid lithium al- 
lowed to detect experimentally the presence of the re- 
laxation scenario predicted by theory on this prototypi- 
cal alkali metal In particular, the data have been 
analyzed within a generalized hydrodynamic approach 
pO| , p^ and, following a memory function formalism [Q, 
it has been possible to detect the presence of two relax- 
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ation processes (additional to the thermal process result- 
ing by the coupling of density and temperature variables) 
affecting the dynamics of the systems in the mesoscopic 
wavevector region JTst . This results, obtained with a phe- 
nomenological ansatz for the memory function originally 
proposed in a teorethical work by Levesque et al. Q , has 
been then deeply discussed pointing out the physical ori- 
gin of such dynamical processes Q . 

Although alkali metals are commonly referred to as 
the "paradigm of simple liquids" , other liquid metals al- 
ways exhibit the characteristic structural and dynamical 
features that can be interpreted within the elementary 
approaches typically utilized when dealing with Lennard- 
Jones or alkali interatomic potentials [ |l5|Jlq| . A remark- 
able example is provided by liquid aluminum. Several 
theoretical and numerical works have been re- 
ported on this system. In fact, the basic common features 
of simple liquids, such as the presence of collective exci- 
tations in the coherent dynamic structure factor below 
Qm and the positive dispersion of the sound speed as- 
sociated to them have been extensively studied in this 
liquid. On the other hand, in this Q region no exper- 
imental data have been reported to provide a real test 
and/or comparison of collective dynamics. Indeed, in 
liquid Al, the high value of the isothermal sound speed 
Co 4700 m/s prevents for kinematical reasons a study 
of the collective dynamics by inelastic neutron scatter- 
ing (INS): below Q,n the available energy transfer is too 
limited to investigate the acoustic excitations properties. 
Accurate INS data have instead been reported at higher 
wavevectors where, however, the single particle response 
is mainly probed (see for example Ref. p9|]). 

In this work we present the study of the coherent 
dynamic structure factor, S{Q,uj), of liquid aluminum 
(T = 1000 K) by Inelastic X-ray Scattering, in a wavevec- 
tor region, 0.05Qm < Q < O.SQ™, of crucial importance 
as far as the dynamical features of the collective motion 
are concerned. 

II. THE EXPERIMENT 

Our IXS experiment has been carried out at the high 
resolution beamline ID 16 of the European Synchrotron 
Radiation Facility (Grenoble, F). The incident beam is 
obtained by a back- scattering monochromator operat- 
ing at the [hhh) silicon reflections, with /i=9,ll. The 
scattered photons are collected by spherically bended Si 
crystal analyzers working at the same [hhh) reflections. 
The total energy resolution function, measured from the 
elastic scattering by a Plexiglas sample, has a full width 
half maximum of 3 meV for /i=9 and of 1.5 meV for 
h=ll: this higher resolution conflguration has been used 
at Q = 1 and 2 nm^^, where the width of the spectral 
features starts to be comparable to the resolution width 
at /i=9. The wavevector transferred in the scattering 



process, Q=2kiSin{6s/2) (with ki the wavevector of the 
incident photon and 6s the scattering angle) is selected 
between 1 nm^^ and 14 nm~^ by rotating a 7 m long 
analyzer arm in the horizontal scattering plane. The to- 
tal (FWHM) Q resolution has been set to 0.4 nm"^ A 
flve analyzers bench was used to collect simultaneously 
flve different Q values, determined by a constant angu- 
lar offset of 1.5° between neighboring analyzers. Energy 
scans have been performed by varying the temperature 
of the monochromator with respect to that of the ana- 
lyzer crystals. Each scan took about 180 min, and each 
spectrum at a given Q was obtained from the average of 
2 to 8 scans, depending on the values of h and of Q. The 
data have been normalized to the intensity of the incident 
beam. The liquid aluminum sample (0.999% purity) has 
been kept in an AI2O3 container with optically polished 
single crystal sapphire windows (0.25 mm thick). The 
sample length, which for an IXS scattering experiment 
is optimal if coincident with the absorption length, has 
been set to 1.0 mm in order to be nearly optimized for 
both h = 9,11 energies (21747 and 17794 eV respec- 
tively). The cell was then host into a molybdenum oven 
in thermal contact with a tantalum foil heated by the dis- 
sipation of about 100 W of power. All the environment 
has been kept in 10~^ mbar dynamic vacuum. 

The IXS spectra I{Q,u}), collected at fixed Q as func- 
tions of the exchanged energy, are reported in Fig. 1. As 
Qm — 25 nm~^, from the raw spectra we can observe 
the relevant region of dispersion of the acoustic mode. 
Initially, the frequency of the inelastic peaks increases 
linearly with Q. Despite an increase with Q of their 
width, the modes are recognizable even at the highest 
wavevectors explored in the experiment. As expected, 
the ratio between the energy loss/gain sides of the spec- 
tra is ruled out by the detailed balance condition. The 
quasi-elastic portion of the spectra, associated to some 
still unrelaxed damping mechanism, shows a linewidth 
that increases with Q, as the characteristic timescale of 
such process suddenly decreases at decreasing the domi- 
nant wavelenght of the observed density fluctuation. 

At the lower Q-values a non-negligible contribution 
coming from the empty container scattering is clearly vis- 
ible (arrows in Fig. 1). This feature basically stems from 
the inelastic scattering of the sapphire windows. As the 
longitudinal phonon velocity in AI2O3 is more than 10000 
m/s, this "spurious" contribution is beyond the energy 
region relevant for liquid Al and does not signiflcantly af- 
fect the main spectral features of interest. However, for 
an accurate quantitative assessment of the spectral shape 
at the lower Q values, in the analysis we have excluded 
the undesired portions around the AI2O3 phonons by cut- 
ting off the two corresponding energy ranges (shadowed 
areas in Fig. 2). 
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FIG. 1. IXS spectra of liquid aluminum at the indicated 
wavevectors. The instrumental resolution is AE = 1.5 meV 
at Q =1 and 2 nm~^ and AE — 3.0 meV elsewhere. The 
total integration time is about 500 s for each point. 



III. DATA ANALYSIS 

The spectra S{Q,oj) have been described within the 
framework of the generalized Langevin equation for the 
density correlator <j)iQ,t) = (pQ{t)p^Q{0)) /(IpqWI^): 

(Q, t) + u;l{Q)(f>{Q, t) + /o Af (Q, t~t')'(t){Q, t')dt' = 
where 

^o(Q) - [kBTlmS{Q)]Q^ - (co(Q)Q)^ (1) 

In these equations, m is the atomic mass, S{Q) the 
static structure factor, cq{Q) the Q-dependent isother- 
mal sound velocity and M{Q^t) the full memory func- 
tion. Recalling that the Fourier transform of (f>{Q,t) 



is S{Q,uj)/S{Q), and introducing the Fourier-Laplace 
transform of the memory function as M{Q,uj) = 
M {Q, uj) + iM'{Q, uj), the previous expression becomes: 



S{Q) 



TT-^uji{Q)M'{Q,uj) 



[LuM'{Q,u;)r 



(2) 



All the details of the microscopic interactions are now 
embodied in the memory function M(Q,t). For the lat- 
ter, we have allowed a two-time relaxation mechanism of 
non-thermal contributions by assuming that 

M(Q,i) = (7-l)c.2(Q)e~^-Q=t (3) 
+ A2(Q) U(Q)e-*/"°('3) + (1 _ ^(Q))e-*/v(Q) 

where r's, and A are the timescales, the total vis- 
cous strength and the relative weight of the two pro- 
cesses respectively. In Eq. (||) the first term comes from 
the coupling between density and thermal fluctuations 
(7 — CpjCy is the specific heats ratio and Dt — n/nCv 
the thermal diffusivity) , while the other two contribu- 
tions describe the relaxation of purely viscous decay 
channels. As is well known, for the latter the familiar 
viscoelastic model assumes a single exponential decay. 
However, as already found in liquid lithium | [r^ , [l^ , this 
simple ansatz is not accurate enough to reproduce the 
details of the IXS spectra in liquid Al. 

Before any discussion of the fitting procedure, we recall 
that the actual scattered intensity is proportional to the 
convolution between the experimental resolution function 
and the true (quantum mechanical) dynamic structure 
factor Sq{Q,uj), affected by the detailed balance condi- 
tion. To account for the latter, we have used the following 
approximation 



(4) 



which connects Sq{Q,oj) with its classical counterpart 
S{Q,u;). 

Summing up, we used Eq.(H,3) as a model function. 
We modified it according to Eq.(4) and, finally, we folded 
it with the experimental resolution. The result has been 
utilized as fitting function F{Q,lu) for the scattered in- 
tensity I{Q, uj): 



F{Q,uj) = E{Q) / R{u-uj')Sq{Q,u:')dw' 



(5) 



where the constant E{Q) depends on each analyzer effi- 
ciency and on the atomic form factor. The parameters 
S{Q) and uj^[Q) are related by the Eq. so that the 
most obvious procedure would be to put the data on an 
absolute scale estimating S{Q) in a fitting-independent 
way and to fix loq{Q) accordingly. Consequently, i) the 
constant E{Q) in Eq. (||) would be E{Q) = 1 ii) the only 
free fitting parameters would be the relaxation times and 
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strengths of the viscous channels (for the thermal quan- 
tities Dt{Q) and ^{Q) it has been utilized their Q ^ 
value). This method has been successfully applied in the 
lithium experiment, where the S{Q) has been estimated 
from the raw spectra utilizing the first two sum rules 
corrected for finite resolution effects . In the present 
case instead, due to the presence of empty cell contribu- 
tions in the tails of the scattered intensity, an estimate 
of the spectral moment would be not reliable. For this 
reason we have chosen to leave E{Q) and Wo(Q) free 
parameters. The best-fitted value of the latter quantity 
has been then used to calculate S{Q), putting the data 
on absolute scale. 



0.60 



0.30 



0.00 
0.40 

Oh 



^ 0.20 

m 0.00 

0.30 



0.15 



0.00 



Q=13.8 nm 




-75 75 

Frequency co (ps^) 



FIG. 2. Selected examples of the fitting procedure. Open 
circle with the error bars are the IXS data, the full lines are the 
best fit to the data. The shadow boxes indicate the energy 
windows that have been cut off' in the fitting session. The 
resolution function is also reported (dotted line). 

Some examples of the quality of the described fitting 
procedure are reported in Fig. 2, while the reliability of 
the normalization method is illustrated in Fig. 3, where 
the insert reports in a log scale the parameter Wg (Q) . The 
value of S{ff) at small Q is seen to be in fair agreement 
with the thermodynamic value (arrow) p8| . 



In Fig. 4 the relative weight A[Q) = l\\{Q) j {l\\(ff) + 
A'^{Q)) is shown. Compared to the hquid lithium, A{Q) 
in liquid Al appears to have a faster decrease with the 
wavevector, indicating a somehow smaller role of slow 
features at finite wavevectors in this system. 
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FIG. 3. Inset: as free fitting parameter. The 

full line is the expected power two behavior. Figure: Val- 
ues of S(ff) deduced by the previously determined a;o(Q) as 
S{Q) = KTQ^/rmjl{Q). The arrow indicates the thermody- 
namic value of 5(0) by compressibility data. 
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FIG. 4. Ratio, A(Q), between the slow and the total relax- 
ation strength. Full dots (•) hthium, open dots (o) aluminum, 
both slightly above the melting point. In order to compare the 
two systems, the exchanged wavevector has been normalized 
to the static structure factor maximum. 

In Fig. 5 we report for both mechanisms the relax- 
ation times and, in the inset, the quantity lji{Q)t{Q), 
where loi{Q) is the maximum of the longitudinal current 
correlation spectra. 
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FIG. 5. Relaxation times and oji{Q)t{Q) values deter- 
mined by the fitting 

The characteristic time t^^ of the microscopic process 
moderately decreases with Q, while the much longer is 
almost constant, except an abrupt initial decrease at very 
low Q which may be possibly due to finite resolution ef- 
fects. From some algebraic manipulations of Eqs. (||,||), 
it is possible to infer the relations between the param- 
eters ruling each mechanism and the main spectral fea- 
tures, i.e. the peak position and the width of the acoustic 
mode, as well as the intensity and the width of the cen- 
tral quasi-elastic contribution. From the inset of Fig. 5 
it can be seen that for the slow process the condition 
^i{Q)Ta{Q) > 1 always holds, so that this mechanism 
contributes mostly to the elastic linewidth rather than 
to the acoustic mode broadening. On the other hand, 
the fast microscopic process is such that uji{Q)t^{Q) < 1 
at low Q, so that in this region it contributes to the acous- 
tic damping. However, at wavevectors Q > 5 nm~^, one 
reaches the crossover uji{Q)t^{Q) w 1, so that the fast 
mechanism drives an effective increase of the sound speed 
toward a limiting, solid-like, instantaneous response. As 
the weight of the fast process exceeds the one of the slow 
one (namely A{Q) « 1 in Eq. (^), cf. Fig. 4), the in- 



crease of the sound speed basically stems only from the 
fast microscopic mechanism, which above the crossover 
threshold is also responsible for the broader tails of the 
quasi-elastic portion of the spectra (tq >> r^). 

In Fig. 6 we report the deconvoluted spectra, i.e. the 
S{Q,uj) as obtained from the fitting procedure, together 
with the relative current correlation spectra J{Q,uj) — 
Lu^ /Q^S{Q,uj). The latter quantity plays an important 
role in the description of the dynamical features of a 
systems because its maxima are related to c; = oJi/Q, 
namely the apparent sound speed of the system. 
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FIG. 6. Deconvoluted spectra of the density correlation 
function S{Q,u!) (left axis) and of the current correlation 
function J{Q,u)) (right axis) as obtained from the fit. 

Finally, in Fig. 7, we compare the Q-dependent 
sound speed in liquid Al (calculated as the maximum 
of the fitting-deconvoluted longitudinal current correla- 
tion function J{Q, ui) = /Q'^S{Q, uj)) with the one pre- 
viously determined in liquid lithium. Suitable reduced 
units (velocities normalized to their respective isother- 
mal values, wavevectors measured in units Qm) have been 
adopted. 
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FIG. 7. Effective sound speed (dots) of aluminum 
(full symbols) vs lithium (open symbols). Data have been 
scaled by Qm -the first sharp diffraction peak and co(0) 
-the isothermal speed of sound in the hydrodynamic limit 
(cq* = 4450m/s;c^' = 4700m/s). co(Q) is also reported 
(— V — ) together with the Coo{Q) values determined by the 
fitting (— A —) and by the pair distribution function and in- 
teratomic potential (lines). The lines connecting the symbols 
are guidelines for the eyes only. 

Both systems are seen to exhibit a clear positive dis- 
persion effect, i.e. an increase of the sound speed at in- 
creasing the exchanged wavevector. Such dispersion, in 
ordinary liquids, proceeds between the adiabatic {Q 0) 
limit, Cs = \/icq, and the unrelaxed, instantaneous value 

Coo((5) = V7t^o(<3) + A2(Q)/Q, where is the total 
strength of the non-thermal processes. In metallic liq- 
uids, instead, due to the high thermal conductivity, the 
condition lui{Q) < DtQ^ holds for Q ^ O.lnm"^, i.e., 
in the IXS Q-window (Q > lnm~^), no region exists 
where the three relaxation processes thermal, a, fj,, are 
simultaneously unrelaxed. As a consequence, the lower 
and upper value of the sound speed are expected to be 
the isothermal co{ Q) (defined in Eq.([l|)) and the partially 

umclaxed c'^ = ^uj^iQ) + AKQ) + AliQ)/Q (20 1. 

To test the reliability of our model, we finally compared 
c^^iQ) with the theoretical value c^^(O) deduced from 
the structure and interatomic potential data ||l^,^. In 
both systems the fitted values of c^'^{Q 0) slightly 
exceed the data for c^{0) reported in literature. A pos- 
sible explanation of this inconsistency may be ascribed to 
the arbitrary choice of the memory function shape, that 
in principle can have more complicated features than the 
double exponential ansatz of Eq. (^). The major draw- 
back of such an assumption is, indeed, the cusp ai t — 0. 
It is reasonable to think that an exponential decay forced 
to represent a more complicated time dependence can 
give an accurate estimate as far as the r is concerned, 
while at short times the lack of a zero second derivative 
inescapably leads to an overestimate of Ml{Q, t — 0), the 
positive dispersion amplitude. A similar effect, i.e. the 



overestimation of the Coo(Q) deduced by the fit has been 
also observed in liquid lithium, where the same memory 
function has been adopted p^. 



IV. CONCLUSION 

In this work inelastic X-ray scattering technique has 
been utilized to study in detail the main features of the 
microscopic dynamics in liquid alumimim in the meso- 
scopic momentum region (Q « 1 — 15 nm^^), i. e. the 
region were the collective properties are dominant. Alu- 
minum is a system where the investigation of the dynam- 
ics using the well-established inelastic neutron scattering 
technique is not possible. The present IXS experiments 
allowed for the determination of the pattern of relax- 
ation processes entering in the density-density memory 
function, and therefore affecting the spectral shape of 
the dynamic structure factor. Specifically, the high qual- 
ity of the data allows to identify univocally the presence 
of three distinct relaxation processes. The first one - 
the usual thermal relaxation- is associated with the cou- 
pling between density and energy fluctuations. As in 
alkali metals, although its relevance is rather small, this 
process cannot be neglected; more interesting, it plays 
a role which is different with respect to the ordinary 
non-conducting liquids. Much more important is the un- 
ambiguous evidence of two well separated timescales in 
the decay of that part of the memory function which 
id associated to the generalized viscosity. In particu- 
lar, the present experiment proves that the traditional 
description of the generalized viscosity by a single time- 
scale (viscoelastic model) cannot account for the detailed 
shape of S{Q,lu), and that it only provides a qualitative 
description of the microscopic dynamics in liquid metals. 

The presence of a two time-scales decay of the memory 
function poses a question of crucial importance: namely, 
their physical origin and, most important, of the fast one. 

The slow process, making use of the terminology used 
to describe glass forming systems, is related to the a- 
relaxation which, in systems capable to sustain strong 
supercooling, is responsible for the liquid-glass transi- 
tion. From a more canonical point of view, this process 
can be framed within kinetic theory in terms of "corre- 
lated collisions" : in mode coupling approaches, the onset 
of these correlation effects is traced back to the coupling 
to slowly relaxing dynamical variables, and specifically 
in the liquid region, to long lasting density fluctuations. 
The quantitative description of this "slow" timescale re- 
quires a full evaluation of the mode-coupling contribu- 
tion at different wavevectors which is beyond th aim of 
the present work. 

As far as the fast process is concerned, the situation 
is more confused and its interpretation is still matter of 
debate. This is particularly annoying because -as shown 
in this paper for liquid aluminum and as demonstrated in 
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the case of liquid lithium- the fast process is indeed more 
relevant than the slow one: it largely controls both the 
sound velocity dispersion and attenuation in the meso- 
scopic Q region. Historically, within a generalized kinetic 
theory approach, the fast initial decay of the density fluc- 
tuation memory function, is traced back to collisional 
events, which are fast, short range and, more important, 
uncorrelated among each others. Within this scheme, the 
short timescale turns out to be associated with the 
duration of a rapid structural rearrangement occurring 
over a spatial range = 2'K/Qm- Although this description 
of the fast process in term of uncorrelated interparticle 
collision is a possible way to qualitative account for the 
dynamical features at short times in liquids and dense 
systems, it seems to be unable to account for similar re- 
sults obtained in glasses pi 22 24|, where the collective 



aspect of the dynamics can not be neglected. 

A possible, and different, approach to describe the ini- 
tial fast decay of the memory function in dense fluids 
and liquids relies on the normal modes ("instantaneous" 
in normal liquids) analysis of the atomic dynamics [p5| , 
an approach that works correctly in the limiting case 
of "harmonic glasses". In this case one uses a frame- 
work (dynamical matrix, etc.) formally similar to the 
one customarily adopted for harmonic crystals; however, 
owing to the lack of translational symmetry of the sys- 
tem, it turns out that the eigenstates cannot anymore 
be represented by plane waves (PW) even at relatively 
small wavevectors p6[ . A scattering experiment, where 
Q is fixed, is equivalent -via the fluctuation-dissipation 
relation- to a response experiment, where one study the 
time evolution of the system after a sinusoidal pertur- 
bation (with period 2ttIQ) is applied to the system it- 
self. In the case of measuring the S{Q,uj), the applied 
perturbation is a density fluctuation. As the PW are 
not eigenstate of the disordered systems, the initial per- 
turbation can be projected along different eigenmodes, 
each one evolving in time with its own frequency. As a 
consequence of this frequency spread the energy initially 
stored in the PW relax towards other PW's with differ- 
ent Q values. This process is exactly what one expect 
in presence of a relaxation process, and its characteristic 
time is determined by the projection of the system eige- 
modes on the PW basis (Fourier transform of the eigen- 
vectors). We deal, therefore, with a mechanism whose ul- 
timate origin is the topological disorder, and not a truly 
dynamical event. Ordinary liquids (such as the one con- 
sidered here) are certainly "disordered systems" and at 
the high frequency considered here {uJTa >> 1), they can 
be considered as "frozen"; therefore the description used 
for glasses can be applied as well. Obviously, the "har- 
monic approximation" and the residual effect of the finite 
value of ujT prevent the " quantitative" application of the 
previously described formalism to the case of normal liq- 
uids. However we expect that -on a qualitative ground- 
the fast relaxation process in liquids can be described in 



term of " high frequency vibrations" (that are a strongly 
correlated atomic motion) , and this point of view appears 
to be in contrast with the (uncorrelated) binary collision 
based description. 

At the present stage the experimental data cannot sup- 
port either interpretation. However, the possibility that 
the phenomenology reported here for a simple liquid can 
be understood within the same mental framework devel- 
oped for more complex liquids and glasses, is certainly 
appealing. 



[1] 
[2] 
[3] 

[4; 

[6: 
[9: 

[lo; 

[11 



[12; 

[is: 

[14 

[is: 
[16: 
[17: 
[is: 

[19 
[20 



J. R. D. Copley and M. Rowe, Phys. Rev. A, 9, 1656 
(1974). 

T. Bodensteiner, Chr. Morkel, W. Glaser and B. Dorner, 
Phys. Rev. A, 45, 5709 (1992). 

C. Morkel and W. Glaser, Phys. Rev. A 33, 3383 (1986); 
A. Stangl, C. Morkel, U. Balucani and A. Torcini, J. 
Non-Cryst. Solids 205-207, 402 (1996). 

P. Verkerk, P. H. K. De Jong, M. Arai, S. M. Bennington, 
W. S. Howells and A. D. Taylor, Physica B 180 & 181, 
834 (1992). 

A. G. Novikov, V. V. Savostin, A. L. Shimkevich, R. M. 
Yulmetyev and T. R. Yulmetyev, Physica B 228, 312 
(1996). 

P. Chieux, J. Dupuy-Philon, J. F. Jal and J. B. Suck, J. 

Non Cryst. Sohds 205-207, 370 (1996). 

H. Mori, Prog. Theor. Phys., 33, 423, (1965). 

D. Levesque, J. Verlet and J. Kurkijarvi, Phys. Rev. A 
7, 1690 (1973). 

J. Casas, D. J. Gonzales and L. E. Gonzales, Phys. Rev. 
B 60, 10094 (1999); J. Casas, D. J. Gonzales, L. E. Gon- 
zales and M. Silbert, preprint of LAM 10, Dortmund 
(1998). 

U. Balucani and M. Zoppi, Dynamics of the Liquid State 
(Clarendon Press, Oxford, 1994). 

T. Scopigno, U. Balucani, A. Cunsolo, G. Ruocco, F. 



Sette, R. Verb eni, Europhys. Lett. 50, 189 (2000); cond 



mat/9911021 



J. P. Boon and S. Yip, Molecular Hydrodynamics 
(McGraw-HiU, New York, 1980). 

T. Scopigno, U. Bal ucani, G. Ruocco, F . Sette, accepted 
by Phys. Rev. Lett. |cond-mat/0009145| 



T. Scopigno, U. Balucani, G. Ruocco and F. Sette, Jour- 



nal of Physics 12, 8009 (2000). cond-mat/0001190 



M.M.G. Alemany, J. Casas, C. Rey, L.E. Gonzales, and 
L.J. Gallego, Phys. Rev. E 56, 6818 (1997). 
W. Gudowski, M. Dzugutov, and K.E. Larsson, Phys. 
Rev. E 47, 1693 (1993). 

T. Gaskell, J. Phys. F 16, 381 (1986) and Phys. Chem. 
Liq. 17, 11 (1987). 

I. Ebbsjo, T. Kinnel and I. Waller, J. Phys. C: Solid State 
Physics, 13, 1865 (1980). 

O.J. Eder, B. Kunsch, J.B. Suck and M. Suda, Journal 
de Physique, 41, C8-226 (1980). 

By a numerical point of view, the difference between 



7 



Coo{Q) and c'^XQ) is of the order of (3%), namely, 
within the statistical errors of the fitting procedure. It is 
worth to point out how the thermal process, contrarily 
to what happens in ordinary liquids, gives also a non- 
viscous damping Sto = {-y — 1) c%{Q)/Dt, responsible for 
an additional Brillouin linowidtli particularly important 
at the lower Q's of our experiment. 
[21] T.Scopigno et al., in preparation. 

[22] G. Ruocco, F. Sette, R. Di Leonardo, G. Monaco, M. 

Sampoli, T. Scopigno and G. Viliani, Phys. Rev. Lett. 

84, 5788 (2000). 
[23] L M. de Schepper and E. G. D. Cohen, J. Stat. Phys. 27, 

223 (1982). 

[24] W. Gotze and M.R. Mayr, Phys. Rev. E 61, 587 (2000). 
[25] G. Seeley and T. Keyes, J. Chem. Phys. 91, 5581 (1989). 
[26] V. Mazzacurati, G. Ruocco and M. Sampoli, Europhys. 
Lett. 34, 681 (1996). 



8 



